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I. INTRODUCTION 



The long range behavior of the Coulomb force causes technical problems in the scattering for 
more than two particles. For instance the 3-body Faddeev kernel develops singularities, which 
deny a direct numerical approach. A way out has been searched in the past by starting with 
a screened Coulomb potential, which for instance in the context of the 3-body problem leads to 
a screened 2-body Coulomb t-matrix. In the limit of an infinite screening radius it is claimed 
in the literature 

EL El El 

that the on-shell 2-body t-matrix approaches the physical one except 
for an infinitely oscillating phase factor, known analytically. Thus removing that factor, called 
renormalization, the physical result can be obtained. 

As a basis for that approach work by Gorshkov Ford 0,0] and Taylor is most 

often quoted. Gorshkov 0, regards potential scattering on a Yukawa potential in the limit of 
the screening radius going to infinity. He works directly in 3 dimensions avoiding a partial wave 
decomposition. He sums up the perturbation series to infinite order. As a result he finds the 
limit for the wave function of a Yukawa potential for an infinite screening radius. That limit 
function equals the standard Coulomb wave function multiplied by an infinitely oscillating phase 
factor. Contrary to what is quoted in Chen he has not achieved the wave function for a Yukawa 
potential at an arbitrary screening radius but only its limiting form. 

The work by Ford 0, relies on a partial wave decomposition. This leads to a very difficult 
technical task to handle the situation, when the orbital angular momentum I is about pR, where 
p is the asymptotic wave number and R the screening radius. This task is left unsolved and the 
infinite sum in I is carried out without controlling the /-dependence of certain correction terms 
depending on R. In other words the correction terms for given I are assumed to remain valid also 
for the 3-dimensional objects. This leaves at least doubts about the rigorousness of that approach. 
The same is true for the investigations of Taylor , where again a partial wave decomposition 
is the basis and the infinite sum over I is carried through without control of its validity for the 
correction terms. 

In such a situation we felt that a rigorous analytical approach for a sharply cut off Coulomb 
potential carried through directly in 3 dimensions is in order. This paper delivers an analytical 
solution for an arbitrary cut-off radius. Further we also provide an exact expression for the cor- 
responding scattering amplitude (equivalent to the on-shell t-matrix). The paper is organized as 
follows. In section [H] the wave function is derived. In section [TTTl the scattering amplitude and its 
limit for vanishing screening is given. These purely analytical results are confirmed by numerical 
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studies presented in section IIVI In the Appendix we regard the much simpler case for s-wave 
scattering. We summarize in section IVl 



II. THE WAVE FUNCTION FOR A SHARPLY CUT-OFF COULOMB POTENTIAL 

Let us regard two equally charged particles with mass m. Then the 2-body Schrodinger equation 
reads 

(_V 2 -p 2 + ^)^(+)(f) = . (1) 

It is well known that in parabolic coordinates 

u = r — z (2) 
v = r + z (3) 

4> = tan~ 1 ^ (4) 

x 

the partial differential equation factorizes and yields the solution 

vj/(+)(j^) = cons t e 1 ^' 7 \Fi(— irj, 1, i(jpr — p ■ r)) (5) 

2 

with Somerfeld parameter r\ = 

Now we switch to a sharply screened Coulomb potential 

e 2 

V(r) = e(R-r)— (6) 
and rewrite ([1]) into the form of the Lippmann-Schwinger equation 



(2tt) 3/2 4^ 



r — r 



This defines uniquely the wave function \£ 1 (r) for a given cut-off radius R. Acting on ([7|) with 
(-V 2 -p 2 ) and using the well known property of the free Greens function in the integral kernel 
one obtains the Schrodinger equation 

2 

(-V 2 - p 2 )^\r) = -m@(R - r )—^\r) . (8) 
Thus for r < R one has to have 

V { +\r) = Ae^Fx {-irj, 1, i{pr - p- r)) (9) 
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with some to be determined constant A. The idea is therefore, to insert that knowledge into the 
Lippmann-Schwinger equation ([7|) leading to 



4 +) ( f ) = — 

R (2vr) 3/2 

- — / d 3 r'——e{R-r')—Ae vp - r iFU-irj, 1, i{pr' - p ■ f ')) . (10) 
Air J \r — r \ r 

If we choose r < R then also the left hand side is known and one obtains the following identity 
Ae l P' r \Fi(— irj, 1, i(pr — p ■ r)) 



(2vr) 3/2 

d 3 r' % —. Q(R-r')^Aj^\F 1 (-ir 1 ,l,i(pr' -p-r')) . (11) 
|r — t* | r 

This provides the factor A. If j4 is known one can determine the scattering amplitude /r defined 
for r — > oo by 

R { ) {2v f 2 

+ — d 3 r'e- ipr - r @(R-r')^-e ip r M-irj, 1, i(pr' - p • f ')) 

r 47r 7 ' 

(2vr) 3/2 r 

It is not difficult using properties of the confluent hypergeometric function to show that the 
corresponding LS equation, for instance for a s-wave, is identically fulfilled as it should. Doing that 
one can read off the corresponding analytical expression for A. That calculation is deferred to the 
Appendix [XI 

The 3-dimensional case is much harder. Let us choose p = z and work with the parabolic 
coordinates. Then dill) turns into 



Ae^-^xFA-in^ipu) = \ rF :S (v - u) 

(2vr) 3/2 

2 p2R f -2R-u' /-27T p ip\r~r'\ 

+ A- \ du'e-^'xFxi-i^ipv!) dv'S v '{—^) d$ —— . (13) 

2 Jo Jo 47r Jo |r-P| 

Since we want to determine just one factor A one value of u and v is sufficient and we choose the 
simplest case u = v = 0. Then the <j>' integration is trivial and one obtains 

1 . e 2 m f 2R , , „ , „ f 2R ~ u ' , , w 1 



A = ~ A ^T / du'iFxi-i^ 1, ipu') / dv'e**" , (14) 

(2vr) 3/2 2 J Q Jo u' + v' 

where we used iFi(—ir), 1,0) = 1. Substituting v! = 2Rx,v' = 2Ry and defining A = A — \^ one 

(2tt) 

obtains 

rl rl-x i 

A = l-ArjT dxxF x (-i<ri,l,iTx) / dye iT?/ (15) 

io Jo x + y 



with T = 2pR. 

Introducing z = iT let us define 

P 1 /* 1 — X "1 

F(z) = l + ^ dxxFA-i^l.zx) dye zy . (16) 

' Jo Jo X + V 

Substituting zx = T,zy = t' we get 

F(z) = l-ir ] [ dnFti-ir],!^) [ 'dr'/-^ . (17) 



o Jo 



Then it follows 

dF(z) irj 
dz z 



e z / dT 1 F 1 (-i V ,l,T)e- T (18) 



o 



dz 2 



' hJr(t) -- % -\e z f ^dr 1 F 1 (-ir / ,l,r)e- T - ^e z [' dr^-ir,, 1, r) e - T 

- — / dTiFi(-ir],l,T)e- z 
z Jo 

1 dF(z) dF(z) if] „ , . „ . 



Consequently 



We add ir]F(z) on both sides 

d 2 F(z) dF(z) 

+ ir]F{z) = iv{F{z) ~ ii?i(_ir? ' 1; z)) • (21) 

The left side put to zero is the defining differential equation for \F\{—ir], 1, z). Thus (f2Tj) is fulfilled 
for 

F(z) = 1 F 1 (-i V ,l,z) (22) 

which also fixes the normalisation. 

A rather lengthy sequence of analytical steps (not given) using an integral representation, re- 
currence relations and further properties of the confluent hypergeometric function yields the same 
result 

Thus we obtain based on (fT5l) 



A = l-A{ 1 F 1 {-i v ,l,iT)-l) . (23) 

The cancellation of A on the left against A on the right is a verification that the LS equation ([1 
at r = is fulfilled, as it should and we end up with the exact relation 



This is valid for any T = 2pR and therefore 

= (2^/^(4, Mr) 6 *^'-*"- * r " * ^ < 25 > 

is exactly fulfilled for r < R, inside the range of the potential. To the best of our knowledge this 
is the first time that this has been achieved. 

At the same time it provides due to (|12p the exact expression for the scattering amplitude fn 
or the on-shell t-matrix element for a sharply cut-off Coulomb potential. This will be dealt with 
in the next section. 



III. THE SCATTERING AMPLITUDE 



The starting point due to (fT2|) with fx = ^ 2 J)3/^ f R 1S 

h = M~) I^Ve-^'^e^'FUl^'-fr')). (26) 

We use the general integral representation of F(a, (3, z) 

1 F x (a,P,z) = C(a,P) j dte zt t a ~ l {l - if~ a ~ x (27) 

where the path T encircles the logarithmic cut between t = and t = 1 in the positive sense and 
the prefactor is 

c < a -fl = r M rff-a)i-^-) - (28) 

Inserting ([27]) into ([26]) yields 

f R = A{-^-)C(-iri,l) J dt^—^) ir) j J d \> e -ipr-?'^jP-? e i(3>r'-p-?')t _ (29) 

The r ' integral is straightforward and one obtains 

f R = A(-^)C(-ir,,l) [ dt{^r^— 



2p 2 a ; v h ; J r t ' t(l-t) 
(1 + e iRt (it sin RVt* + 2(l ^a _ & _ 
V V + 2(1 - t)a V ^ > >> 



(30) 



where a contains the dependence on the scattering angle 9 



q 

a = 1 — p ■ f = 2sin 2 — , (31) 

C(—ir), 1) = ^e nv , and R = pR. The "1" in the bracket does not contribute since 

cft(— r M 1 , = . (32) 
r v t ' t(l-t) K J 



Thus we obtain the intermediate result 

a* m ^^f^-^in 1 jmsinR^ + 2(1 -t) 



f R = -A^-C(-i V ,l)[i [ dt(\±) 
ap Jr t 



t l-t ^/t 2 + 2(1 - t)a 

dt( lzlyv^_ e iRt cos ^ t 2 + 2(1 _ t)a] _ (33) 
p t i(l — t) 

In the following we choose the path of integration T as depicted in Fig. [T] with small circles 
around t = 1 and t = of vanishingly small radius e and two straight integration lines between 
t = e and t = 1 — e above and below the logarithmic cut. The phases are : arg(t) = and 
ar<7(l — i) = 7r for £ = 1 + e. The rest follows by continuity: arg(l — t) = 2ir along the upper rim of 
the cut, arg(t) = 2tt along the lower rim and arg(l — t) = 3tt back again at t = 1 + e. The phase 
of — p does not change after a full sweep of T, of course. 

In this manner the integrals in (j33H can be split into 4 pieces. Let us define 



Then 



B = / dtfl^ly* 1 i& r - ^\A 2 + 2(l-*)<* _ cosR^/t 2 + 2(1 -t)a 
1 { t ] l-t { ^2 + 2(1 _ t )a t 



B = [ + f \ [ + f • (35) 

J zero Je J one Jl—e 



It simply follows 



1-e 



(H + 1 1 C cft = (1 - e" 2 ^) / e <ft(I_*)*J_L_ e i& 



,.sm J R v / t 2 + 2(l-t)a cosR^/t 2 + 2(1 - t)a 



(36) 



y/t 2 + 2(1 - t)a * 
In order to remove the pole singularities at t = and t = 1 we split the integration interval into 
two parts 

l-e 1-1/2 pl-e 

dt= dt+ dt . (37) 

Je Jl/2 

Of course the value 1/2 could be replaced by any number a between t = e and t = 1 — e without 
changing the result. 
Thus 



''"dt = i / 1/2 ^(l-tr-^V- Sm ^ + 2(l-t)a 



yft 2 + 2(1 -t) 



o 



+ i r^dtQ . r -U^^nR^t 2 + 2(1 -t)a 
Jl/2 ^t 2 + 2(1 -t)a 

- J dtt-^-^l - t) ir i- l e iRt cosR^t 2 + 2(1 - t)a 



rl-e . 

- / dt(l - ty r >-H- ir i- 1 e iRt co8Ry/t 2 + 2(1 - t)a . (38) 

A/2 



r-l-e 
1/2 

Now we perform partial integrations such that e — ► can be taken: 

-1-e Z" 1 / 2 • - • ikt sinR^f t 2 + 2{l- t)a 



A Jo y / t 2 + 2(l- t)q 

l iV y/t 2 + 2(1 - t)a 1/2 



+ 1 r e dm f^d ( ^ iVri fy sinRy/t 2 + 2{l-t)a , 
"7 A/2 dt K y/t 2 + 2(1 -t)c 



_ r J_ r iTjM _ tY^e^cosRWt 2 + 2(1 - t)a|V 2 

— 17] 

1 /" 1/2 , d 

' ' )n) 

*»7 Jo 
1 



+ _ / dtt~ iri ^((l - tp- l e iRt cosRJt 2 + 2(1 - t)c 
vq J dt 



[-L-(l - t )^r^- l e iRt cosR^t 2 + 2(1 - t)a\)j* 



— 17] 

+ — f dt(l - tt' — ir^e^cosRJt 2 + 2(1 - t)a)] . (39) 
it] J 1/2 dt 



After some lengthy algebra one obtains 



r- dt = ■-,-.^,■4°^ 

A "7 "7 ?7 Vl/4 + a 

Jo 7 Jo 



+ « 

/o J ir l Jo 

- f 1 dt(l-tyn-^e iRtS -^(t-a)-i I' dtil-trr^e™?^ 
VJi/2 T Ji/2 7 



+ dt(l - tyn- iT >- 2 e iAt cosR-f 

11] Jl/2 

i> r 1/2 - 1 fl/2 _ / — o 

+ -[- / dirtl - tyi-^cosR-y + - / dit^l - tyi-^sinRr/ 

?7 Jo i Jo 7 

+ i f 1 d t{i-tyH-^ s -^L + C dtii-tyn-^^L^ 

A/2 7 A/2 7 7 

- / cft(l - i)^^" V^cosi^ + - / dt(l - f) , 't- i '- 1 e^s!nfl7— ] 

A/2 « Jl/2 7 



(40) 



with 7 = y/t 2 + 2{\-t)a. 

It is straightforward to evaluate the two integrals around t = and t = 1: 

1-t ^ 1 iAt smgvgTgl^ga _ cosR^t 2 + 2(1 - t)q 
zero t 1 — i y 7 * 2 + 2(1 — tja t 

-icosRV2^e- ir >^(l - e~ 2 ^) , (41) 

J one t ' l-t y/t 2 + 2(1 -t)a i 
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I e *7(i _ e -an,) . (42) 

7? 

The e-dependent terms cancel in (HDD (multiplied by (1 — e _27r7? ) ), (|4ip and (|42|) as they should 
and one obtains the finite result 

B = (1 --^^^VV^ + a | . f 1/ \ t{1 _ t yv-i t -iv e im s J^LL 
V a/1/4 + a J 7 

+ ^— *: [ 1/2 dtt-^(l - ty^e^cosR-f - - f 1 dt(l-m-^e ikt ^^(t-a) 
w) Jo VJi/2 T 

f 1 dt(l - t yn t ~^-i e im s j!^L + fH+1 f' dt h _ ifit^^cosBq 

Jl/2 7 ^ Jl/2 



















r 






1/2 




R 




+ 




[" 




V 










+ 


l J 


r 






1/2 






1 



/■1/2 - 1 fl/2 - / _ „ 

/ diT^l - t) ir '- 1 e im cosR 1 + - / dtr^l - tP^e^sinR-y 
Jo 1 Jo 

+ f dtil-tfH-^e 

Jl/2 



7 

iv iRt sinRy f 1 j iv iAt cosRn/t--a 



7 Ji/2 7 7 

/ dt(l - tfH-^e^cosRy + - [ dt(l - t)^"*^ 1 e^sinflT— ]] . (43) 

Jl/2 * Jl/2 7 

This together with ()33|) - ()36| ) is an exact expression for the scattering amplitude for an arbitrary 
cut-off radius R. 

But of course we are interested only in its asymptotic limit R — > oo. 

It is advisable to introduce e± = e %R ^^ and to rearrange (|43p . We regard first the pieces 
explicitly proportional to R in (f4"3j) 



2r] J t 1-t J t 1-t 7 

+ dt(l^r e+ ~ e - + C dt^ir e -±- +e - t - a 



1/2 * ' 7 Ji/2 t 7 7 



dt(—)^-(e + + e_) - / df(-_)«>-( e+ - e_) ] . (44) 

1/2 * * Jl/2 * * 7 

Leading terms will arise from the boundaries of integration t = 0, t = 1/2, and t = 1, where 
the t = 1/2 contributions have to cancel in the total expression. We use the standard method of 
steepest descent [9J and expand around boundaries of integration. For example at t = 



e± = e ±tiiv^ e iitt(lTV§)(l + 0(t)) (45) 
and corresponding expressions for the remaining parts of the integrand. One obtains 



2?7 L Jo < 1 - 1 Jo * 1 - i 7 

J_(_ ie 5»Jr(l - ^) J R^(e 2 ^ n l (1 - sffl^" + g-^fcmf^ + sin ^)<i) . ( 46 ) 
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Correspondingly we proceed at the upper limit of integration t = 1 and it turns out that the e+ 
part decreases as O(-i) and only the e_ part survives as 



2rj J t 7 J t 7 7 



t ' i v ' / v £ ' i v ^ 7 

-(2i?)- i,? (sin 2 ^)- ir 'e^r(l + i? ? ) . (47) 

The remaining pieces resulting from the integration limits t = 1/2 yield 

o /-1/2 rl/2 f — n 

£[- / - t) iT >-\e + + e_) - / tttt-^Cl - t) < "- 1 (e+ - e_)— - 

2?? 7 / 7 

, e + — e_ /" , j j e+ + e_ t — a 



A/2 7 /l/2 



'1/2 7 Ji/2 7 7 

- / dt(l-trt- iri - 1 e iRt (e + + e^)- [ dt(l - t)^t~ ir '- 1 (e + - e_)^^] 

/l/2 A/2 7 

-» --— — e^sinRy/TJI+a ■ (48) 

This cancels exactly against the first term in ([43 p after multiplication by (1 — e~ 27rr? ), as it should. 

The terms in (|43p not directly proportional to R decrease like O(i). Finally the contributions 
from the interior of the integration intervals decay faster as can be seen by deforming the path of 
integration into the upper half plane, where e^- is exponentially damped. 

Thus we are left with the leading asymptotic expression 

B -» (l-e~ 2 ^)— ef"( -T(l-iri)R ir i( e 2iksin i {I - sin-)^ + e ~ 2iksin i (I + sin-) ir > ) 
2r\ 2 2 

+ 2{2R)- ir '(sin 2 ^)-^T(l + i V ) ) . (49) 



This is now to be combined with (|33p . Using (|28l) . 

ee e 2lr7 ° (50) 



r(l -irj) 
and the asymptotic form of A 

e-^T-^T(l + ir]) (51) 

based on the asymptotic form [lo| 

±iira — a z a—f} i 

i%,A^ L , +^-^ + Q(^) , (52) 



r(/3-a) r(a) >| 



we get 



/ fl = -(2fl)" A " e 



ot r7 (sin tt) 



2p sin 2 1 
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Now the physical Coulomb scattering amplitude is 



AcW = -^^|pe 2 - (54) 
2p svn l \ 



and we end up with 



f R = (2R)~ 2ir 'A c (e)+ V fl (e 2 ^"f( (1 f!!^ )^ + e -2iw|( (i±!!2g 



Apsin 2 1 2 2 



_ |- e -2ir??ra(2_R) _ ^i-qlnsin 2 1— 2i<7o ^iRsin^+irjln 1 "^"^ _|_ e -2iRsin^+irjln 1+s ^ w 2 ^ (55) 

The first term is the result expected from the literature and references therein. As has 
shown, the diverging phase factor e~ 2l ® R ^ in case of an often used form of screening the Coulomb 
potential 



V R (r) = ^e~^ n (56) 

n 

using the prescription of [lj turns out to be 

$ fl (p) = V [ln(2 P R) - C/n] (57) 

with the Euler number C. For n — > oo one recovers the sharp cut-off, which we consider in this 
paper. This expectation for the screening limit agrees with the first term in (|55|) but not with the 
necessity of adding a second term. Therefore the derivations in the literature based on partial wave 
decomposition must be incomplete. Whether this is also true for a finite value n in (|56p remains 
to be seen. 



IV. NUMERICAL RESULTS 



We performed a number of numerical tests to check the basic points in the derivation of the 
sharp cut off Coulomb wave function (|25p and the asymptotic scattering amplitude (|55p. 

First we checked numerically how well the solution (|24|) fulfills equation (|15|) . In Table [J the 
left and right sides of (|15p are shown for a number of cut-off radii R for pp scattering with E l ^ h = 
13 MeV. The right side was obtained by a direct two-dimensional numerical integration over x and 
y. The very good agreement up to four significant digits is seen. 

We also compared at the same energy the exact expression for A as given in (|24p with its 
asymptotic form (|5ip at a number of screening radii. The results are shown in Fig. [2] and Table ITTl 
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The oscillating behavior seen in real and imaginary parts of exact A (solid lines in Fig. [2]) gradually 
diminishes with increasing cut-off radius R. These oscillations are absent in the asymptotic form 
for A (dashed lines in Fig. [5]). The asymptotic form for A approaches its exact value at R ~ 50 fm 
as can be seen in Fig. [2] and in the third column of Table HT1 where the ratio of A/A approx is given. 

To check the quality of our renormalization factor (|55p we applied it directly to the numerical 
solutions of the Lippmann-Schwinger equation for the sharp cut off Coulomb potential with different 
cut-off radii. 

In the case of a short-ranged potential V two-body scattering is described by the solution of 
the Lippmann-Schwinger equation 

T{z)=V + V—^—T{z), (58) 
Z — Hq 

where V is the two-body potential, Hq is the free Hamiltonian and T(z) the transition operator. 
In momentum space Eq. ()58[) takes the form of an integral equation for the matrix elements of the 
transition operator (q \ T(z) \ q) = T(q',q). In this equation matrix elements of the potential 
V are used (q' \ V \ q) = V(q',q). In our case both V(q',q) and T(q',q) depend only on the 
magnitudes q' =\ q' \, q =\ q\ and the cosine of the angle between q and q ', q ' ■ q: 

V(q',q) = V(q',q,q' • q) (59) 

T(q',q)=T(q',q,q' ■ q). (60) 
(Note we dropped the dependence on the parameter z.) As a consequence the Lippmann-Schwinger 



equation can be written as a two-dimensional integral equation 

oo 1 
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T(q',q,x') = ^v(q' ,q,x',l) + J dq"q" 2 J dx"v(q', q",x',x") ^ 7T T(g // ,g, x"), (61) 



o 

where 

/ / / \ 

V 



(q',q,x',x) = J d(pV(q',q,x'x + \/l — x' 2 \/l — x 2 cos ip) (62) 
Q 

and m is the reduced mass of the system. 

For the sharply screened Coulomb potential of the range R considered in this paper 

rtf,«,«) = ^ 1 -'% Wfl) . (63) 

where Q = yq 1 2 + q 2 — 2q' qx. However, the integral over ip in Eq. (|62p cannot be carried out 
analytically. 
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It is clear that V(q',q,x) shows a highly oscillatory behavior, especially for large R. Thus 
solving the two-dimensional equation (|61|) is a difficult numerical problem. We were interested in 
solutions for positive energies where 

z = E c . m . + ie = ^ + ie. (64) 
m 

We solved (|6ip by generating the corresponding Neumann series and summing it up by Pade which 
is a very reliable and accurate method. Usually six iterations were fully sufficient. In each iteration 
the Cauchy singularity was split into a principal-value integral (treated by subtraction) and a 5- 
function piece. We used 120 or 140 (/-points and 150 or 190 x-points. The (/-integral points are 
chosen in the definite interval (0, q), where typically q= 50 fin . In order to obtain directly the 
on-shell t-matrix element T(qo, qo,x, E c _ m _) we added q = qo to the set of (/-points. To better control 
the behavior of the transition matrix element for small scattering angles also x = 1 was added to 
the set of appoints. A typical run required less than 9 minutes on 256 nodes (1024 processors) on 
the IBM Blue Gene/P parallel computer at the Jiilich Supercomputing Centre. 

In Figs. [3][5] we show with dash-dotted line the real and imaginary parts of the transition 
amplitudes Ac(0) = —2ir 2 mT(qo,qo,cos9) for sharp cut off Coulomb potential pp scattering at 
E l p ab = 13 MeV and a number of cut-off radii R = 10 and 20 fm (Fig. [3]), R = 40 and 80 fm 
(Fig. Hj), and R = 100 and 120 fm (Fig. Ej). With increasing cut-off radius a development of strong 
oscillations in the scattering angle dependence for the real parts of the numerical solutions is clearly 
seen. These oscillations follow on average the real part of the pure Coulomb amplitude given by 
(|54p and shown by the solid line. The imaginary parts of the numerical solutions are totally off from 
the imaginary part of the pure Coulomb amplitude and have even an opposite sign. Now applying 
to the numerical solutions the asymptotic renormalization factor from (|55jl dramatically improves 
the agreement (dotted lines in Figs. OS]). Not only the oscillations in the real parts are practically 
removed and the pure Coulomb and renormalized amplitudes are practically overlapping but the 
renormalization brings also imaginary parts into agreement with the exception of very forward 
angles. When one desists to use the asymptotic expansion for /r and instead calculates it exactly 
according to (|33p . (|34|) and (|43p than the ratio ^j^y provides the exact renormalization factor. 
Performing exact renormalization of the numerical solutions provides very good agreement between 
imaginary parts of the numerical and pure Coulomb amplitudes also at the very forward angles 
(dashed line in Figs. [3]E])- 

We also checked how important are the two additional terms in the renormalization factor of 
(I55p . To this aim we renormalized the numerical solutions with the standard form of the renor- 
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malization factor, given by the first term in (|55[) . In Fig. [6] solid (red) lines show the amplitude 
renormalized in this way. It is clearly seen, that restricting to the standard form of the renormal- 
ization factor it is not possible to reach the physical amplitude. Standard renormalization reduces 
slightly oscillations in the real part of the numerical solution and changing the sign of the imagi- 
nary part invokes in it large oscillations. So after standard renormalization strong oscillations are 
present both in the real and imaginary parts and fails totally. 

V. SUMMARY 

The renormalization method for a screened on-shell Coulomb t-matrix enjoys a widespread 
use; see for instance [l! As pointed out in the introduction the underlying mathematical 
considerations leave room for doubts. To shed light on that issue we regarded potential scattering 
on a sharply cut-off Coulomb potential directly in 3 dimensions, avoiding obstacles in the infinite 
sum of angular momenta. The idea was to use the Lippmann-Schwinger equation which uniquely 
defines the wave function including its boundary conditions. Inside the range of the potential it is 
the standard Coulomb wave function multiplied by an unknown normalisation factor. Using that 
form also on the left side of the Lippmann-Schwinger equation for radii smaller than the cut-off 
radius determines that normalisation factor uniquely. Based on that we succeeded analytically 
to determine the normalisation factor and thus obtained in this manner the exact analytic result 
for the wave function. This also allowed us to derive the analytical expression for the scattering 
amplitude in the limit of infinite cut-off radius. The connection to the standard Coulomb scattering 
amplitude A c (6) turned out, however, to be different from the standard form used widely in the 
literature and is given in (|55p . Our form consists of two terms, one of which is the standard 
one, e~ 2ir]ln2pr A c (6) . To that, however, is added a new expression which is singular at 6 = and 
9 = tt. These analytical results are fully backed up by accompanying numerical investigations. 
Our renormalization factor brings in a very good agreement between the strongly deviating and 
oscillating numerical solution of the Lippmann-Schwinger equation with the sharp cut off Coulomb 
potential and the exact Coulomb amplitude. The standard renormalization factor fails completely. 
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APPENDIX A: S-WAVE POTENTIAL SCATTERING FOR A SHARPLY CUT-OFF 
COULOMB POTENTIAL 

The (reduced) wave function for s-wave scattering obeys the Lippmann-Schwinger equation 

m f R , , , , e 2 



0(+)(r) =sin(>r) / dr'e ipr > sm(p r< ) — ^ + \r') (Al) 

P Jo r 

with i"<(>) the smaller ( greater) of r,r'. Inside the potential range (f)( + \r) has to have the form 

</>(+) (r) = AF (pr) (A2) 
where F$(pr) is proportional to the standard Coulomb wave function 

F (pr) = pre ipr F(l + in, 2, -2ipr) . (A3) 
Inserting (|A2|) into (|A1[) yields 

[ r 1 I ^ i 1 

(j)( + \r) = sin(pr) — 2r]A(e ipr / dr' sin(pr') — F^ipr') + sin(pr) / dr'e ipr — F^ipr')) 

Jo r J r r 

= sin(W) - ^A( e vr f dr'e 2ipr 'F(l + in, 2, -2ipr') 
R 2^ 

- e~ ipr [ dr'e 2ipr 'F(l + irj,2,-2ipr')-e ipr I dr'F(l + in,2,-2ipr')) . (A4) 

One faces two types of integrals, which can be solved using the following properties of the 
confluent hyper geometric function: 

1 d 

F(l + in, 2, -2ipr) = 1, -2ipr) (A5) 

2prj dr 

F(l + in, 2, -2ipr) = / dte- 2iprt (— — ) ir ? (A6) 

2ttt] J r 1-t 

with the path T given in section IIIIl and 

F(irj, 1, -2ipr) - F(l + irj, 1, -2ipr) = 2iprF{\ + in, 2, -2ipr) . (A7) 

One obtains 

r i 

/ dr'Ffl + in, 2, -2ipr') = (Ftin, 1, -2ipr) - 1) (A8) 

Jo 2 CT 

/ dr'e 2ipr 'F(l + it;, 2, -2»pr') = — (1 - e 2ipr F(l + in, 1, -2ipr) . (A9) 

Jo 2r ?P 
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Therefore 

0(+)(r) = sinfpr) - ^—(-^- S m(pr)e 2ipR F(l + in, 1, -2ipR) - -^e ipr F(l + in, 2, -2ipr)) 

2z 2r?p ' 2i]p 

= sin(pr)(l - Ae 2ipR F(l + in, 1, -2ipR)) + ,4pre ipr F(l + ir/, 2, -2ipr) 

= <p( + \r) + sin(pr)(l - Ae 2ipR F{l + ir/, 1, -2ipR)) . (A10) 



Consequently the LS equation (|A1|) is identically fulfilled, as it should and one obtains an explicit 
condition for the constant A: 

1 - Ae 2ipR F{l + ir), 1, -2*pi?) = (All) 

or 

„-2ipR 

A = — . (A12) 

F(l + irj,l,-2ipR) v ' 

Inserting this result into (|A2|) the exact s-wave function for a sharply cut-off Coulomb is obtained 

-2ipR 

* ( ,(r) = FiHWP) F '*' f(1 + 2> <A13) 



It obeys the LS equation (jATj). 

The asymptotic behavior r —* oo, which provides the scattering phase shift 5r(p), is given 
through the LS equation and we read off from (|A4|) 



0(+) ( r ) ^ sin (p r ) _ gipr^/ ( Al4 ) 

with 

pR 

A' = 2rjpA / dr' sm(pr')e ipr ' F(l + ir], 2, -2ipr) . (A15) 
•/ o 

At the same time this yields 

e 2i5 R (p) = 1 _ 2iA , _ ( Alg ) 



Using (|A8[) and (|A9|) again gives 

A' = — (e 2ipR F(l + ir/, 1, -2ipR) - Ffin, 1, -2ipR)) (A17) 
2z 

and consequently 

e 2i5 fl (p) = j_ _ A(e 2ipR F(l + in, 1, -2ipi?) - F(ir/, 1, -2ipi?)) . (A18) 
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The interest lies now in the limit R — > oo. We use (|A12p and the asymptotic form (|52p of F 
and obtain 

e 2i8 R (p) _^ e 2ia -2iri\n(2pr) (A19) 

or 

5 R (p) ^ a - vH^Pr) ■ (A20) 

Of course this result is well known and can be trivially obtained by matching the interior 
Coulomb wave function to the free one containing 5r(p). 

We performed this exercise to explicitly demonstrate that the LS equation (|Aip is indeed iden- 
tically fulfilled for arbitrary r below the cut-off radius R. In the 3-dimensional case we succeeded 
analytically to do this only for the special value r = 0, though it is valid for any r < R, and were 
forced to verify the general case numerically. 
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TABLE I: The left and right sides of at E l p ab = 13 MeV (77 = 0.0439, p = 0.3959 fm" 1 ) and different 
screening radii R. 



R [fm] 


A 1 


- AtjT J 1 dxxFii-ir}, 1, iTx) £ x dy^^- 

1 j \ 111 / j \j " x~\-y 


0.5 


(0.98301, -0.00166) 


(0.98301, -0.00166) 


1 


(0.96724, -0.00634) 


(0.96724, -0.00634) 


5 


(0.91933, -0.08246) 


(0.91933, -0.08246) 


10 


(0.92770, -0.10294) 


(0.92770, -0.10294) 


20 


(0.91961, -0.13491) 


(0.91961, -0.13491) 


50 


(0.91606, -0.17185) 


(0.91606, -0.17185) 


100 


(0.90960, -0.20061) 


(0.90960, -0.20061) 


500 


(0.89376, -0.26439) 


(0.89377, -0.26439) 


1000 


(0.88528, -0.29140) 


(0.88528, -0.29140) 


5000 


(0.86250, -0.35307) 


(0.86252, -0.35307) 



TABLE II: The exact value of A as in (|24[) (left column), asymptotic form given by (f5Tj) (middle column) 



and their ratio (right column) at EJ 1 — 13 MeV for different screening radii R. 











R [fm] 


A 


A 

si-approx 


A I ' A a pp rox 



0.1 (0.99653,-0.00007) 

1.0 (0.96724,-0.00634) 

2.0 (0.94156,-0.02249) 

3.0 (0.92536,-0.04353) 

5.0 (0.91933,-0.08246) 

10.0 (0.92770,-0.10294) 

20.0 (0.91961,-0.13491) 

50.0 (0.91606,-0.17185) 

100.0 (0.90960,-0.20061) 



(0.92852, 0.07999) (1.06534,-0.09185) 

(0.93186,-0.01402) (1.03784, 0.00881) 

(0.93100,-0.04233) (1.01035, 0.02178) 

(0.93010,-0.05888) (0.99389, 0.01611) 

(0.92855,-0.07970) (0.99040,-0.00380) 

(0.92570,-0.10788) (1.00152, 0.00551) 

(0.92199,-0.13596) (0.99731, 0.00074) 

(0.91578,-0.17289) (1.00009 0.00115) 

(0.91011,-0.20064) (0.99946,-0.00008) 
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FIG. 2: (Color online) The real (top) and imaginary (bottom) part of A as a function of the screening radius 
R at E l p ab — 13 MeV. The solid (black) line represents the exact expression given in (|24|) and the dashed 
(red) line shows the asymptotic form as given in (|51[) . 
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FIG. 3: (Color online) The real (top) and imaginary (bottom) part of Ac(9) = —2ir 2 mT(qQ,qo,cos9) as 
a function of cos 9 for R— 10 fm (left panel) and 20 fm (right panel) at E l p ab = 13 MeV. The dash-dotted 
line represents a direct numerical prediction (without any renormalization) . The dotted line shows Ac{&) 
with inclusion of the asymptotic renormalization factor given in (|55|) and the dashed (red) line is for Ac (9) 
with inclusion of the exact renormalization factor obtained from (f3"3")) . (j3"4"|) and (|4"3")l (see text). The solid 
line represents the pure Coulomb amplitude given in (|54|) . Note that the dashed, dotted and solid lines 
practically overlapp with exception of very forward angles for imaginary part. 
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FIG. 4: (Color online) The same as in Fig. [3] but for R= 40 fm (left panel) and 80 fm (right panel). 
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FIG. 5: (Color online) The same as in Fig.[3]but for R= 100 fm (left panel) and 120 fm (right panel). 
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FIG. 6: (Color online) The real (top) and imaginary (bottom) part of Ac(9) = —2ir 2 mT(qQ,qo,cos9) as 
a function of costf for R= 20 fm (left panel) and 100 fm (right panel) at E l p ab = 13 MeV. The dotted line 
represents a direct numerical prediction (without any renormalization) . The solid (red) line shows Ac{&) 
with renormalization factor e - 2,i ;M 2 pfi) an d the dashed line represents the pure Coulomb amplitude given 
in HI. 
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